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Abstract 

We analytically and numerically investigate the long-term, i.e. averaged over one full 
revolution, orbital effects of the non-isotropic percent mass loss m/m experienced 
by several transiting hot Jupiters whose atmospheres are hit by severe radiations 
flows coming from their close parent stars. The semi- major axis a, the argument 
of pericenter u and the mean anomaly M experience net variations, while the 
eccentricity e, the inclination / and the longitude of the ascending node remain 
unchanged , on average. In particular, a increases independently of e and of the 
speed V esc of the ejected mass. By assuming \m\ < 10 17 kg yr , corresponding to 
\rh/m\ < 10~ 10 yr" 1 for a Jupiter-like planet, it turns out d ~ 2.5 m yr -1 for orbits 
with a = 0.05 au. Such an effect may play a role in the dynamical history of the 
hot Jupiters, especially in connection with the still unresolved issue of the arrest 
of the planetary inward migrations after a distance a > 0.01 au is reached. The 
retrograde pericenter variation depends, instead, on e and V csc . It may, in principle, 
act as a source of systematic uncertainty in some proposed measurements of the 
general relativistic pericenter precession; however, it turns out to be smaller than 
it by several orders of magnitude. 
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1 Introduction 



As shown first bv l 2 1 iMescerskiil (118971 ). the equation of motion for a body 
acquiring or ejecting mass due to interactions with its environment i 



dv 

~dt 



F 



(1) 



with respect to some inertial frame K. In it, rh 
the external forces, and 



dm/dt, F is the sum of all 



u 



(2) 

is the velocity of the escaping mass with respect to the barycenter of the body: 
v is the velocity o{3 the center of mass of the body with respect to K, whild_] 
v esc is the velocity of the escaping particle with respect to K. If the mass loss is 
isotropic with respect to the body's barycenter, then the total contribution^! 
(rh/m)u vanishes. Notice that eq. ([Q) fulfils the Galilean invariance under a 
transformation v — > v = t> — V to another inertial frame K moving with 
constant velocity V with respect to K. 

An interesting, real physical scenario in which eq. (JT|) is applicable is rep- 
resented by those transiting exoplanets whose atmospheres are escaping! be- 
cause of the severe levels o f energetic radiations, coming from their very close 
parent stars, hitting them ( lEhrenreich and Desertl . 1201 ll ). For those few tran- 
siting planets which are observed in the ultraviolet it is possible to gain in- 
formation about the size and mass-loss rate of their evaporating upper atmo- 
spheres. It turns out t hat typical values of the ey aporation rate for HD 209458b 



( ICharbonneau et all. |2000| 



HD 189733b (Bouchv et al. 



Henry et al.l . |2000[ ). unoff icially named "Os iris' 



2005h. and WASP-12b flHebb et all l2009h 



in the range ( lEhrenreich et al.l . l2008t iFossati et al.l . l2010t iLinsky et al.l . 12010 



are 



2 See also Sommerfeld ( 1952 Y_ Razbitnava ( 1985 ). Hadjidemetriou ( 19631 ) and the 
review by lHadjidemetrioul (jl967l ). 

3 It is just the case to remind that an erroneous form of eq. (fTl), in which — v appear s 



instead of u, was not rarely used in literature. See, e.g., Plastino and Muzziol ( 1992 ) 
for a discussion of such a topic. Suffice it to say that it would be non- invariant under 
Galilean transformations. 

4 More precisely, it is the velocity of that point of the body which, at each moment, 
coincides with the body's center of mass. It does not include the geometric shift of 
the center of mass caused by the mass loss. 

5 To avoid confusions, it may be useful to remark that in literature v esc is sometimes 
denoted by u, while no specific notations are used for the relative velocity in eq. 
©. Cfr., e.g.. lRazbitnaval (Il985l h 

6 Strictly speaking, it is to be intended as (^ i rhiUi)/m, where the label i refers 
to the ith point of the surface from which mass is escaping. 

7 See |http://www.spacetelescope.org/news/heic0403/| on the WEB. 
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Lecavelier Pes Etangs et all boiol ; Ividal-Madiar et all 12003 . 12004 120081 1 

3 x 10 14 - 10 15 kg yr" 1 = 10 10 - 10 11 g s -1 ; (3) 



m 



however, lEhrenreich and Desertl (120111 ) estimate the mass-loss rates of all de- 
tected transiting planets finding an upper bound as large as 



\m\ < 3 x 10 17 kg yr^ 1 = 10 13 g s~ x . 



(4) 



Several, still unresolved issues are connected with giant exoplanets like the 
aforementioned ones. Indeed, they all orbit at less than 0.1 au from their 
parent stars: thus, they certainly could not have formed there, so that an 



inward migration from more distant locations should have occurred (ILin et al. 



19961 ; iLubow and Ida! 120101 ) . Why did such hot Jup iters stop migrating at such 
distances? What was the fate of those gaseous giants that migrated further 
in, if any? Is the evaporation responsible of the apparent lacking of exoplanets 
closer than 0.01 au to their host stars? 

Studying the orbital consequences of eq. ([1]) may help to shed some light 
on such important open problems. In Section [2] we work out, both analytically 
and numerically, the long-term variations of all the standard six Keplerian 
orbital elements induced by eq. ([T]). Then, we apply the results obtained to 
some specific exoplanetary scenarios. Section [3] summarizes our findings. 



2 Orbital effects 



In the case of an evaporating transiting exoplanet orbiting a close Sun-like 
star, F is the usual Newtonian gravitational monopole. Let us write down eq. 
([!]) for both the planet p, with mass m . and its hosting s tar s, with mass M, 
with respect to some inertial frame K (IRazbitnayal . 119851 ) : 



dt 



dv s 
dt 



n \ m / 



M ' \M 



U 



(5) 



where 



and 



GMm 



GMm 



r 3 V^P 5 



r 3 \ r s fpj , 



r = \r p - r s \ . 



(6) 



(7) 
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In eq. (JSj) we took into account the possibility that also the star experiences 
a mass loss due to internal ph ysical processes. In the cas e of the Sun, it is 
estimated to be of the order of (jSchroder and Smith! . 120081 ) 



M 

— = -9 x 10- 

M 



14 yr" 1 ; 



about 80% of such a mass-loss is due to the core nuclear burning, while the 
remaining 20% is due to average solar wind. In any case, it can be considered 



(see, e.g., C 


Jvlden ( 


L884); Stromeren (1 


903) 


Jeansl {'. 


L924): 


Armellinil (] 


935 


): 
); 


Vescan 


(19371. ] 


939h: Jeansl f 


1961 


): 


Hadiidemetriou 


(1963) 


Omaro^ 


/ (1 


964 


Hadiidemetriou 


(1966) 


Kevorkian and Cole 


1966) 


Kholshevnikov and Fracassini 


(196£ 


0; 


Verhulst and Eckhaus 


(1970) 


: Kurvshev and Perov ( 


L981); 


Depri 


;J1983); 


Polvakhova 


(1989): 


Prieto and Docobo 


(1997a b): 


Li et al. 


(2003 


); ] 


( 


20081, 


2009 


;] 


-lahoma et al. ( 


2009): 


!orio 


(2 


OlOaf)) is devoted on treating the motion 



in a two-body system with various kinds of time-dependent isotropic mass loss 
affecting the Newtonian monopole itself; for a recent treatment of such a topic 
with the same method of the prese nt paper and a discussion of some of the 



approaches present in literature, see lloriol (j2010al ). Concerning the evaporating 



planet, let us assume for it a constant percent decrement fn/m of its mass; it is 
especially true over timescales of the order of one orbital period P\>, if the orbit 
is assumed almost circulai0- For a typical Jovian mass m ~ mj = 1.899 x 10 27 
kg and the figure of eq. (T4J) we have 



\m/m\ < 1.7 x 10~ 10 yr 



(9) 



Moreover, its mass decrement is clearly non-isotropic with respect to its center 
of mass, so that (rh/m)u p does not vanish. As a result, from eq. ([5]) and 
eq. (P) it can be o btained the equation for the relative planet-star motion^ 
( Razbitnaya . 19851 ) 

dv /i frn\ 

— ^r+(— u p , 10 

mj 



dt 



where we defined 



[i = G(M 



m) 



r p - r s , v 



Moreover, u p can be conveniently expressed as the difference between the 
escape velocity V esc with respect to the star and the velocity v of the planet 



It should be not th e case ( Iro and Deming . 20ld : Ehrenreich and Desertl . 201 ll ) 



for, e.g., HD 80606b dNaef et al 



20011 ) whose orbit has an eccentricity as large 
as e = 0.93. Recall that it is a numerical parameter determining the shape of the 
Keplerian ellipse: < e < 1, with e = co rresponding t o a ci rcle. 



^ Cfr. with the classification in Table I of iB^bitnavJ S): eq. (TTO]) falls within 
the A.II.b.9 (Seeliger) or B.II.b.20 (Fermi) cases depending on the relative sizes of 
m and M. 
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with respect to the star. Clearly, V esc is radially directed from the parent star 
to the planet, i.e. 



^"csc ^esc-Rj R 



(12) 



anc0 (IVidal-Madiar et all l2004h 



V PSC ~ 10 4 m s" 1 . 



(13) 



In order to evaluate the magnitude of u p , let us note that in most of the 
considered exoplanets the orbits are almost circular, so that it can be posed 



+ VL -2(vV e 



V 2 



(14) 



by assuming that a Jupiter-sized planet is at a distance^ 1 "! a of the order of 
0.05 au from its Sun-like parent star, it turns out 



v rsj Jtl „ 10 5 m s _1 , 
V a 



(15) 



so that 



Up ~ v ~ 10 5 m s . 



(16) 



Such figures imply that the mass-escaping term in eq. (ITU]) can be considered 
as a small perturbation with respect to the Newtonian monopole A^; 
indeed, from eq. (jHJ) and eq. (flBT) turns out 



Am ~ 2 m s" 



(17) 



Arn ~ 5 x 10 -13 m s 



Thus, it can be treated with the standard perturbative techniques like, e.g., the 
Gauss equations for the variati on of the osculating Keplerian orbital elements 
(IBrouwer and Clemencd . Il96ll ) which are valid for any ki nd of disturbing ac- 
celeration A, irrespectively of its physical origin. See also lOmarovl (119621 ) for 
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Actually, such a figure comes from the escape speed from the planet q e 



^j2Gm/R p : in order to have V^ sc one should also add the radial component of the 
planet's motion which, however, is negligible because of the assumed low eccentricity 
of its orbit (see below, eq. ([20]) ). 

11 Here a is the semi-major axis of the Keplerian ellipse: it defines its size. 
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their use in the variable mass case. The Gauss equations are 



da 
dt 



nVT- 



eA R sin / + A T 



de _ VT- 
dt m 



{A R sinf + A T [cos/+i(l-j|)]} 



dJ_ 

dt 



na\f\ 



\r^ A N (£) cosn, 



da 

dt 



. } -r — 7 An ( l ) sin 



d-uj 
dt 



sfl- 



■A fl co8/ + A r (l + 5Jsin/ 



f2sin 2 (|)f, 



dM 
dt 



n-^A 

na 



R 



+ cos/ 



r dQ 
dt 



n = y/i/a 3 is the unperturbed Keplerian mean motion, e is the 
p = a(l — e 2 ) is the semi-latus rectum, / i g 12 1 the inclination 



In eq. (ED 

eccentricity, 

of the orbital plane to the reference {x,y} plane chosen, Q is the longitude 
of the ascending node, u is the argument of the pericenter, M. is the mean 
anomaly! 13 1. / is the true anomaly, u = u + / is the argument of latitude, 



and Ar, At, A^ are the radial, transverse and out-of-plane components of the 
disturbing acceleration A, respectively. To this aim, can be written 



Am — I 



(V csc -v R )R-v T f 



(19) 



where v r ,vt are the 14 1 radial and transverse components of the planet's ve- 
locity evaluated onto the unperturbed Keplerian ellipse 



Vt 



nae sin E 
1— ecos E ' 



(20) 



iVT- 



l—e cos E ' 



12 The parameters /, f2,w defines the spatial orientation of the Keplerian ellipse, 
which, in the unperturbed case, changes neither its size nor its shape: they can be 
thought as the three Euler angles fixing the orientation of a rigid body in the inertial 
space. 

13 



ii 



It is connected with the time of pericenter passage t p through M. = n(t — t p ). 
Here T denotes the unit vector along the transverse direction, i.e. orthogonal to 



R. 
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In eq. (120]) . E is the eccentric anomaly 15 1. From eq. (fT9l) it is inferred that 
An = 0, being Ar and At the only non-zero components of A^; they are 



.4 



R 



nae sin £* 
1— e cos E I ' 



(21) 



■ e cos E ' 



The long-term effects of eq. (|T9|) can be straightforwardly worked out af- 
ter an integration of the right-hand-sides of the Gauss equations in eq. (|T8|) . 
evaluated onto the unperturbed Keplerian ellipse where ix is to be assumed as 
constant, over one orbital revolution by means of 



dt 



e cos E 



n 



dE: 



(22) 



we assume rh/m constant over one orbital period. If the mass variation oc- 
curred at fast rates with respect to the orbital frequency, it should explicitly 
be treated as a specific function of time in eq. ffTUl) and eq. ffT^l) . They are 



da 
dt 



dQ 
dt 



dt 



"2« (S) 



— = 

dt u ' 



dt 



dM _ _ fm\ 
dt \ m J 



rn\ V csc Vl~e 2 



3V CS 



(23) 



Notice that eq. ( 123]) are mathematically exact in the sense that no simplifying 
assumptions concerning e were adopted in the calculation. On the other hand, 
from a physical point of view it is necessary to require moderate values for 
the eccentricity since, otherwise, it would not be possible to consider rh/m 



15 It can be regarded as a parametrization of the usual polar angle 9 in the orbital 
plane. 
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constant over the integration performed over one full orbital revolution. It 
is worthwhile noticing that, even in the case of almost circular orbits, the 
long-term variations of eq. (I23p can be considered as secular changes only 
over timescales in which it is possible to assume fn/m as constant; otherwise, 
slow time-dependent modulations, depending on the exact law of variation of 
m = m(t), occur. Here we just recall that the long-term variations caused by an 
is otropic mass loss of the parent star were worked out with the same approach 



in 



Ioriol ( 2010a ): it turned out that all the osculating Keplerian orbital elements 



remain unchanged, apart from the osculating semi-major axis and eccentricity 
which undergo long-term variations 



«iso — 2 ( 1-e ) ( M 



■ \ |§ ) ( i + e ) • 



(24) 



This outwardly counter- intuitive result, and the lack of actual contradiction 



with the occurrin g expansion I of the true orbit for a mass decrease, are 



JXDc 

fully discussed in lloriol ( l2010al ) with abu ndance of exp licative pictures. The 



expressions of eq. (124"]) were also used by lloriol (l2010bl ) in the framework of 
dark matter studies in our solar system. In any case, we notice that, given the 
order of magnitude of isotropic mass losses in typical Sun-like main sequence 
stars (cfr. with eq. (jSJ)), such effects are smaller than those investigated here 
by some orders of magnitude. 

From eq. (|23|) it turns out that the semi-major axis a undergoes a temporal 
variation proportional to rh/m; it is independent of the eccentricity and of 
V esc . The eccentricity is, instead, left unaffected. As expected, if rh < the 
size of the orbit increases since a > 0. 



We qualitatively checked our analytical results by numerically integrating 
the equations of motion of eq. (flUjl . in cartesian coordinates, for a fictitious hot 
Jupiter experiencing a constant mass loss during its motion along a Sun-like 
parent star, and tracing the osculating Keplerian ellipses at two consecutive 
pericenter passages. The results are displayed in Figure [T]for an almost circular 
initial orbit, and in Figure [2] for a highly elliptical initial orbit; they have the 
unique purpose of effectively displaying the mathematical agreement with the 
analytical results, and the percent mass loss was purposely greatly exaggerated 
just to this aim. Given the merely illustrative purpose of the figures displayed 
here, just aimed to numerically confirm the qualitative features of the effects of 
eq. (123]) . longer time intervals for the numerical integrations are unnecessary. 



16 Indeed, from eq. (27) and eq. (37) of Ioriol ( 2010al ) it can be inferred that the 
secular change of the star-planet distance is dA/dt oc — {^M/M \ a(l — e), so that 
the orbit actually expands for M/M < 0, as expected. 
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-0.06 



-0.04 



-0.02 

x(AU) 



0.02 



0.04 



0.06 



: 10 4 m s 1 , in the field of a Sun-like star. The initial 
conditions are x = a (l - e ), yo = 0,z = 0, x = 0, y = n a J j^, z = 0, with 
0.05 au, eo = 0.01, P\> = 2tt/hq = 4.08 d. The numerical integration is over 
2Pb = 8.16 d. Red dashed line: osculating Keplerian ellipse at the first peri- 



Fig. 1. Black continuous line: numerically integrated trajectory of a fictitious Jo- 
vian-type planet experiencing a purposely exaggerated constant percent mass loss 
rh/m = -0.025 d~\ with V e 

ao(l 

a 
At 

center passage (t = 0). Blue dotted line: osculating Keplerian ellipse at the second 
pericenter passage (t = 4.85 d). It can be clearly noticed that the eccentricities of 
both the osculating Keplerian ellipses remain constant, contrary to their semi-major 
axes which, instead, increase. 

The increase of the semi-major axis and the constancy of the eccentricity 
of the osculating Keplerian ellipses ar e apparent. A nalogous pictures for the 
isotropic mass loss case can be found in lloriol (j2010al ). Also from a quantitative 
point of view the agreement with out analytical results is excellent. Indeed, in 
Figure [3] we depict the change of the semi-major axis over a Keplerian orbital 
period obtained from a numerical integration of the equations of motion of 
a fictitious hot Jupiter in cartesian coordinates with and without a constant 
mass-loss o f^l rh/m = —0.00025 d _1 ; both the integrations use ao = 0.05 au 
and eo = 0.1. The overall shift is just equal to the one which can be inferred 
from eq. (123|) . Figure HI obtained for eo = 0.8, confirms the analytical finding 



In order to make a meaningful comparison with our analytical, perturbative re- 
sults we choose a value for \rh/m\ small enough to make the extra-acceleration of 
eq. (|19p much smaller (5 orders of magnitude) than the Newtonian monopole over 
the interval of the numerical integration. 
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-0.12 -0.1 -0.08 -0.06 -0.04 -0.02 
x(AU) 



Fig. 2. Black continuous line: numerically integrated trajectory of a fictitious Jo- 
vian-type planet experiencing a purposely exaggerated constant percent mass loss 
m/m = —0.025 d -1 , with V csc = 10 4 m s -1 , in the field of a Sun-like star. The initial 

conditions are x = a (l - e ), yo = 0,z = 0, x = 0, y = n a yjj^, z = 0, with 
ao = 0.05 au, eo = 0.8, P^ = 2-k/tiq = 4.08 d. The numerical integration is over 
At = 2Pb = 8.16 d. Red dashed line: osculating Keplerian ellipse at the first peri- 
center passage (t = 0). Blue dotted line: osculating Keplerian ellipse at the second 
pericenter passage (t = 4.53 d). It can be clearly noticed that the eccentricities of 
both the osculating Keplerian ellipses remain constant, contrary to their semi-major 
axes which, instead, increase. 



o.oooi 

0.00008 
0.00006 

a 

a 

0.00004 
0.00002 


12 3 4 

t(d) 

Fig. 3. Difference between the semi-major axes of a fictitious Jovian-type planet 
in the field of a Sun-like star computed from the numerically integrated equa- 
tions of motion of eq. (|10p with and without a constant percent mass loss 
m/m = —0.00025 d _1 and V csc = 10 4 m s" 1 . The initial conditions are 

x = o (l - e ),yo = 0,z = 0,x = 0,y = n a Q Jj^,Zo = 0, with a = 0.05 au, 
eo = 0.1, -Pb = 27r/no = 4.08 d. The numerical integration is over At = P^, = 4.08 
d. The net increase of a is apparent, with an overall magnitude equal to the one 
which can be inferred from d = —2a(m/m) of eq. (|23p obtained perturbatively. 
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Fig. 4. Difference between the semi-major axes of a fictitious Jovian-type planet 
in the field of a Sun-like star computed from the numerically integrated equa- 
tions of motion of eq. (|10p with and without a constant percent mass loss 
rh/m = —0.00025 d _1 and V esc = 10 4 m s _1 . The initial conditions are 

x = o (l - e ),yo = 0,z = 0,x = 0,y = n a J±±^,z = 0, with a = 0.05 au, 
eo = 0.8, Pb = 27r/no = 4.08 d. The numerical integration is over At = = 4.08 
d. The net increase of a is apparent, with an overall magnitude equal to the one 
which can be inferred from a = —2a(rh/m) of eq. (|23p obtained perturbatively. Cfr. 
with Figure El 

of eq. f[2"5j) that the net change in a is, actually, independent of e. In Figure 
we show that also the numerically integrated variation of the eccentricity 
does not exhibit any secular trend, in agreement with eq. ( 123|) . According to 
eq. ( )23|) . the pericenter and the mean anomaly experience secular precessions 
which depend on V esc and on e. Such a behavior is qualitatively confirmed 
by Figure |6] displaying the result of a numerical integration of a fictitious hot 
Jupiter obtained for suitably chosen values oirh/m, V esc , eo; once again, it has 
just illustrative purposes. 

Our result about the semi-major axis is, in principle, important since it 
yields a physical mechanism which counteracts the inward migration of the 
exoplanet after it reaches a distances small enough to trigger an effective mass 
depletion. The typical figures previously obtained yields 

a ~ 2.5 m yr" 1 , (25) 

i.e. an orbit with a = 0.05 au increases its size by 2.5 m after each year if the 
planet following it loses mass at a percent rate of about \mjm\ ~ 10" 10 yr _1 ; 
by postulating that m/m remains almost constant over the aeons, we would 
have 

Aa ~ 1.1 x 10 10 m = 0.075 au (26) 

after a time span At = 4.5 Gyr. Thus, the overall orbital evolution of close 
exoplanets may be affected by their mass depletion. Concerning the retrograde 
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2 

t(d) 



Fig. 5. Difference between the eccentricities of a fictitious Jovian-type planet 
in the field of a Sun-like star computed from the numerically integrated equa- 
tions of motion of eq. (jlOp with and without a constant percent mass loss 
-0.00025 d _1 and V eRC = 10 4 m s _1 . The initial conditions are 

0,z 



0,x 



= 10 4 
0, m = n a 



m/m = 

eo = 0.1, Ph = 2ir/riQ = 4.08 d. The numerical integration is over At = = 4.08 



i^a, z = 0, with a = 0.05 au, 



d. The absence of a net increase of e is apparent, in agreement with e 
(|23l) obtained perturbatively. 



of eq. 



pericenter precession of eq. ( 123]) . it may, in principle, be viewed as a potential 
source of systematic bias in some proposed attempts to measure the general 
relativistic pericenter precession (IMiralda-Escudel. 120021 ; iHeyl and Gladman . 



20071 ; IJordan and Bakosl . 120081 ; iPal and Kocsisl . 120081 ) ; anyway, it turns out to 
be smaller by several orders of magnitude than the relativistic one for the 
typical values of the relevant parameters adopted so far. 



3 Summary and conclusions 



We investigated the orbital consequences of a non-isotropic mass depletion 
affecting the evaporating atmosphere of a typical hot Jupiter in a close (a = 
0.05 au) orbit along a Sun-like parent star. 

Analytical perturbative calculation showed that the semi-major axis under- 
goes a long-term increment which is independent of both the escape velocity 
of the atmosphere and of the eccentricity of the orbit, which, instead, remains 
constant. Such results were confirmed, both qualitatively and quantitatively, 
by numerical integrations of the equations of motion in cartesian coordinates. 
By assuming \rh\ < 10 17 kg yr _1 , corresponding to \rh/m\ < 10 -10 yr _1 for a 
Jovian mass, it turns out that an increase of a few meters per year occurs for 
the semi-major axis. Such a mechanism may be important for the dynamical 
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-0.08 -0.06 -0.04 -0.02 

x(AU) 



Fig. 6. Black continuous line: numerically integrated trajectory of a fictitious Jo- 
vian-type planet experiencing a purposely exaggerated constant percent mass loss 
m/m = —0.0025 d _1 , with V esc = 10 7 m s _1 , in the field of a Sun-like star. The ini- 
tial conditions are x = a (l - e ),y = 0, z = 0, x = 0, y = n a ^j^, z = 0, 
with ao = 0.05 au, eo = 0.8, Pb = 2-k/tio = 4.08 d. The numerical integration is 
over At = 2P^ = 8.16 d. Red dashed line: osculating Keplerian ellipse at the first 
pericenter passage (i = 0). Blue dotted line: osculating Keplerian ellipse at the sec- 
ond pericenter passage (t = 2.56 d). The retrograde precession of the pericenter of 
the osculating Keplerian ellipses is apparent. 



evolution of hot Jupiters, especially in connection with the issue of their in- 
ward migrations which should bring them close to their hosting stars from the 
supposedly remote regions of their formation. Indeed, it may explain, at least 
in principle, why such a phenomenon seems to cease after the planet reaches 
a distance of not less than 0.01 au. It is just the case to notice that such a 
behavior would be enhanced, in principle, by the isotropic mass variation M of 
the hosting star as well, although it may be some orders of magnitude smaller 
for typical Sun-like values of M/M . . 



The pericenter of the orbit of the evaporating planet experiences a long- 
term retrograde precession which, however, is not a concern for some proposed 
detections of the general relativistic pericenter precession since it is several 
orders of magnitude smaller. 
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